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ABSTRACT 

Detached, inspiraling and semi-detached, mass-transferring double white 
dwarf (DWD) binary systems are both expected to be important sources for 
the proposed space-based gravitational-wave detector, LISA. The mass-radius 
relationship of individual white dwarf stars in combination with the constraints 
imposed by Roche geometries permit us to identify population boundaries for 
DWD systems in LISA's "absolute" amplitude-frequency diagram. With five key 
population boundaries in place, we are able to identify four principal population 
sub-domains, including one sub-domain that identifies where progenitors of Type 
la supernovae will reside. Given one full year of uninterrupted operation, LISA 
should be able to measure the rate at which the gravitational-wave frequency / 
and, hence, the orbital period is changing in the highest frequency subpopulation 
of our Galaxy's DWD systems. We provide a formula by which the distance to 
each DWD system in this subpopulation can be determined; in addition, we show 
how the masses of the individual white dwarf stars in mass-transferring systems 
may be calculated. 

Subject headings: accretion, accretion disks - binaries: close — gravitational 
waves — stars: white dwarfs 



1. Introduction 

Double white dwarf (DWD) binaries are considered to be one of the most promising 
sources of gravitational waves for LISA,^ the proposed Laser Interferometer Space Antenna 
(Faller & Bender 1984; Evans, Iben, & Smarr 1987; Bender 1998). If, as has been predicted 
(Iben & Tutukov 1984, 1986), close DWD pairs are the end product of the thermonuclear 



^http: / /lisa. nasa.gov 
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evolution of a sizeable fraction of all binary systems, then DWD binaries must be quite 
common in our Galaxy and the gravitational waves (GW) emitted from these systems may be 
a dominant source of background noise for LISA in its lower frequency band, / ^ 3 x 10^^ Hz 
(Hils et al. 1990; Cornish & Larson 2003). DWD binaries arc also believed to be (one of the 
likely) progenitors of Type la supernovae (Iben & Tutukov 1984; Branch et al. 1995; Tout 
2005) in situations where the accreting white dwarf exceeds the Chandrasekhar mass limit, 
collapses toward nuclear densities, then explodes. Because its instruments will have sufficient 
sensitivity to detect GW radiation from close DWD binaries throughout the volume of our 
Galaxy, LISA will provide us with an unprecedented opportunity to study this important 
tracer of stellar populations and it will provide us with a much better understanding of the 
formation and evolution of close binary systems in general. Clearly, a considerable amount 
of astrophysical insight will be gained from studying the DWD population as a guaranteed 
source for LISA. 

Broadly speaking, our Galaxy's DWD binary population should be dominated by sys- 
tems that are in two distinctly different evolutionary phases: An "inspiral" phase, during 
which both stars are detached from their respective Roche lobes; and a semi-detached, "sta- 
ble mass-transfer" phase during which the less massive star fills its Roche lobe and is slowly 
transferring mass to its more massive companion. While DWD binaries may encounter other 
interesting evolutionary phases - for example, a phase of so-called common envelope evo- 
lution, or a phase of rapid, unstable mass transfer - the inspiral and stable mass-transfer 
phases are expected to dominate the population because they are especially long-lived. It 
should be noted that we already have a handle on the size of the galactic population of 
DWD binaries from optical, UV, and x-ray observations. In the immediate solar neighbor- 
hood, there are 18 systems^ (Nelemans 2005; Anderson 2005; Roclofs 2005) known to be 
undergoing a phase of stable mass transfer (AM CVn being the prototype) and the ESO SN 
la Progenitor SurvcY (SPY) has detected nearly 100 detached DWD systems (Napiwotzki 
et al. 2004b). At present, orbital periods and the component masses for 24 detached DWD 
systems have been determined (see Table 3 of Nelemans et al. (2005) and references therein), 
five of which come from the SPY survey. 

During both of these relatively long-lived evolutionary phases, a DWD system's orbital 
period (and corresponding GW frequency) changes on a time scale that is governed by the 
rate at which angular momentum is being lost from the system due to gravitational radi- 
ation, that is, the so-called "chirp" time scale (see the discussion associated with Eq. 17, 



^Three models (Cropper et al. 1998; Wu et al. 2002; Marsh & Steeghs 2002) have been proposed to 

determine the nature of two controversial candidate systems (RX ,10806+15 and V407 Vul) out of these 18, 
which can change the number of known AM CVn systems between 16 and 18. 



-3- 



below). While the system is detached, the orbital separation slowly decreases so the sys- 
tem should emit a GW signal with a characteristic "chirp" signature, that is, the frequency 
and amplitude of the GW signal should monotonically increase with time. During a phase 
of stable mass transfer, however, the orbital separation steadily increases so the GW sig- 
nal should exhibit an inverse-chirp character where by its frequency and amplitude should 
steadily decrease with time. The primary objective of our present study is to analyze the 
imprint that both of these relatively long-lived phases of evolution will have on our Galaxy's 
DWD binary population, as viewed by LISA. 

LISA'S capabilities as a GW detector arc usually discussed in the context of the log(/?.) — 
log(/) diagram, where the GW frequency / is measured in Hertz, and the GW amplitude 
/i is a dimensionless "strain" (generally quoted per \/Hz, reflecting the frequency resolution 
of the data stream). A useful analogy can be drawn between this amplitude- frequency 
diagram and the astronomy community's familiar color-magnitude (CM) diagram. Directly 
from photometric measurements, astronomers can produce a CM diagram that is based 
on the apparent brightness (the apparent magnitude m) of various sources. However, a 
determination of the intrinsic brightness of each source must await the determination of the 
distance r to each source and the corresponding conversion of each measured value of m to 
an absolute magnitude M. LISA's measurement of log(/i) for a given astrophysical source 
is analogous to a measurement of m; it only tells us how bright the GW source appears to 
be on the sky. A determination of the intrinsic brightness of each GW source must await a 
determination of the source distance r and the corresponding conversion of each measured 
value of (the apparent brightness) /i to a quantity that represents the "absolute" brightness of 
the source. For LISA sources, the relevant quantity (analogous to M) is \og{rh). Astronomers 
realize that the underlying physical properties of stars, their evolution, and their relationship 
to one another in the context of stellar populations can only be ascertained from a CM 
diagram if M, rather than m, is used to quantify stellar magnitudes. By analogy, it should 
be clear that the underlying physical properties of DWD systems, their evolution, and their 
relationship to one another in the context of stellar populations can be ascertained only if the 
observational properties of such systems are displayed in a \og{rh) — log(/) diagram, rather 
than in a plot of log(/i) versus log(/). For this reason, our discussion of DWD systems will 
be presented in the context of this more fundamental, but rather under-utilized, "absolute" 
amplitude-frequency domain. 

As we investigate the evolution of DWD systems across LISA's "absolute" amplitude- 
frequency domain, we will utilize a simplified description of the two long-lived evolutionary 
phases mentioned above. We will assume (1) all orbits are circular; (2) the orbital frequency 
is related to the orbital separation via Kepler's third law; (3) the spin of both stars can be 
ignored so that each system's total angular momentum is given by the point-mass expression 
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for orbital angular momentum; (4) the total mass of each system is conserved; and (5) angular 
momentum is lost from each system only via the radiation of gravitational waves and that 
the rate of angular momentum loss is correctly described by a quadrupole radiation formula. 
Simplification ^^4 means, for example, that after the low-mass white dwarf comes into contact 
with its Roche lobe, we will assume that each DWD system evolves along a "conservative" 
mass transfer (CMT) trajectory, and simplification jj^b means that we will be ignoring effects 
that might arise due to direct- impact accretion (Marsh et al. 2004; Gokhale et al. 2006). A 
more thorough analysis that removes some or all of these simplifications is likely to provide 
additional valuable insight into the evolution of DWD populations; Stroeer et al. (2005), 
for example, have expressed concern that the time-rate-of-change of the GW frequency for 
many of LISA's most interesting sources will not be correctly interpreted without a proper 
treatment of tides. In this context, our analysis should be viewed as an important first step 
in what is likely to be a long-term study of the evolution of the DWD binary population 
across LISA's "absolute" amplitude-frequency domain. At the outset, we acknowledge that 
our understanding of this subject has benefitted significantly from the insight of others whose 
work has preceded ours; most notably, we recognize the insightful publications by Paczyhski 



(1967), Faulkner (1971), Evans, Iben, & Smarr (1987), Webbink & Iben (1987), Marsh et al. 
(2004), and Gokhale et al. (2006). 



In the quadrupole approximation (Peters & Mathews 1963; Thorne 1987; Finn & Cher- 
noff 1993), the time-dependent gravitational- wave strain, h{t), generated by a point mass 
binary system in circular orbit has two polarization states. The plus and cross polarizations 
of h{t) generically take the respective forms, ^ /i+ = ^norm cos[0(i)] and — /inormSin[0(i)], 
where the time-dependent phase angle. 



where 0o is the phase at time t = 0, f = florh/Ti is the frequency of the gravitational wave 
measured in Hz, florh is the angular velocity of the binary orbit given in radians per second. 



^Throughout this paper when we refer to experimental measurements of h, we will assume that the binary 
system is being viewed "face on" so that the measured peak-to-peak amplitude of the two polarization states 
are equal and at their maximum value, given by /inorm- If the orbit is inclined to our line of sight, the 
inclination angle can be determined as long as a measurement is obtained of both polarization states as 
shown, for example, by Finn & Chernoff (1993). Because our discussion focuses on Galactic DWD binaries, 
we ignore the effects of cosmological expansion. 



2. 



Parameterization 




(1) 
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and the characteristic amphtude of the wave, 

_ G_ ml^M,M,a^ 

(Mi + M2) ' ^ ^ 

where G is the gravitational constant, c is the speed of hght, r is the distance to the source. 
Ml and M2 arc the masses of the two stars, and a is the distance between the stars. If the 
principal parameters of the binary system do not change with time, then / and /inorm will 
both be constants and the phase angle (f) will vary only linearly in time, so the source will emit 
(monochromatic) "continuous-wave" radiation. If, however, any of the binary parameters 
— Ml, M2, a, or Qorb — vary with time, then /inorm and/or / will also vary with time 
in accordance with the physical process that causes the variation. Here we will only be 
considering the evolution of DWD systems in which the basic system parameters vary on a 
timescale that is long compared to 1//. 

We know from the mass-radius relationship for white dwarfs (see the discussion associ- 
ated with Eq. 12, below) that the less massive star in a DWD binary will always have the 
larger radius. Therefore, in a DWD system that is undergoing mass transfer, we can be 
certain that the less massive star is the component that is filling its Roche lobe and is trans- 
ferring (donating) mass to its companion (the more massive, accretor). With this in mind, 
throughout the remainder of our discussion we will identify the two stars by the subscripts d 
(for donor) and a (for accretor), rather than by the less descript subscripts 1 and 2, and will 
always recognize that the subscript d identifies the less massive star in the DWD system. 
This notation will be used even during evolutionary phases (such as a GR-driven inspiral 
phase) when the two stars are detached and therefore no mass-transfer is taking place. Fur- 
thermore, we will frequently refer to the total mass of the system, Mtot = + Ma and the 
mass ratio, 

which will necessarily be confined to the range < g < 1 because M^ < Ma- Also, it will 
be understood that the limiting mass for either white dwarf is the Chandrasekhar mass, 
Mch = 1.44Mo. 

As mentioned above, throughout this investigation we will assume that Kepler's third 
law provides a fundamental relationship between the angular velocity and the separation of 
DWD binaries, that is, 

^l.-^. (4) 

a 

Relation (4) allows us to replace either QqA or a in favor of the other parameter in Eq. (2). 
Furthermore, we will find it useful to interchange one or both of these parameters with the 
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binary system's orbital angular momentum 



Joich — 



(^) Q. (5) 



where, 

is the ratio of the system's reduced mass to its total mass. For our future discussion, it is 
useful to express the gravitational wave amplitude /inorm and frequency / in terms of Jorb 
and Q. So, 



^norm — , P ' 



/ = -^f^- (8) 

^ "^orb 

We note as well that the so-called "chirp mass" of a given system (Finn & Chernoff 1993) 
is obtained from Mtot and Q via the relation, 

M = MtotQ'/' . (9) 



3. Evolution of DWD Binciries in the Amplitude- Frequency Domain 

3.1. Trajectories and Termination Boundaries 

Detached DWD binaries slowly inspiral toward one another as they lose orbital angular 
momentum due to gravitational radiation. Assuming that Mtot and q remain constant during 
this phase of evolution, Eqs.(7) and (8) can be combined to give. 



9^Tv'^ /(^A/f , \ 5 1 1/3 

TT ^'fj^^ K'f2^ ^ 5.38 [x5y2]l/3 ^g^g^g^ (iQ) 



C2 

where the dimensionless mass parameter. 



1 + g. 



has been defined such that it acquires a maximum value of unity in the limiting case where 
Md — Ma — Mch', otherwise, < K < 1. (We note that in the limiting case of K = 1, the 
chirp mass of the system is M. — 0.871Mch = 1.25Mq.) From expression (10), we see that 
the trajectory of an inspiraling, detached DWD binary in the amplitude-frequency diagram 
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can be determined without specifying precisely the rate at which angular momentum is lost 
from the system. Specifically, because d\ii{rh-aorm) / d\n f = 2/3, trajectories of inspiraling 
DWD binaries will be straight lines with slope 2/3 in a plot of log(r/?.norm) versus log/. 
By way of illustration, inspiral trajectories for systems having three different total masses 
(2.4 Mq, 1.4 Mq, and 0.8 Mq) but all having the same q = 2/3 mass ratio are displayed in 
the top panel of Figure 1. 

The detached inspiral phase of the evolution of a DWD binary will terminate when the 
binary separation a first becomes small enough that the less massive white dwarf fills its 
Roche lobe. From Egglcton's mass-radius relationship for zcro-tcmpcrature white dwarfs, as 
quoted by Verbunt & Rappaport (1988) and also in Marsh et al. (2004), we know that the 
radius of the donor is. 



^ = 0.0114 

Rq 



f M, X -2/3 / M, X 2/3n ^2 r / M, N "2/3 ^ f M,yh "2/3 



where Mp = 0.00057 Mq. Furthermore, from Eggleton (1983) we find that the Roche-lobe 
radius Rl is, 



Rj 



0.49 q 



2/3 



0.6 g2/3 + in(i + gl/3) 



(13) 



The orbital separation - and the corresponding GW amplitude r/inorm and frequency / - at 
which the inspiral phase terminates can therefore be determined uniquely for a given donor 
mass Md and system mass ratio q by setting R^ — Rl and combining expressions (12) and 
(13) accordingly. The termination points of the three inspiral trajectories — marked by plus 
symbols in the top panel of Figure 1 — have been calculated in this manner. The curve 
connecting the sequence of plus symbols in Figure 1 traces out the locus of points that define 
the termination points of the detached inspiral phase of numerous other DWD systems that 
have mass ratios q — 2/3 but that have values of M^ot ranging from 2.4 Mq to 0.06 Mq. 

As a DWD system fills its Roche lobe and starts transferring mass to its companion, 
it evolves to lower frequencies and amplitudes. Without knowing the precise rate at which 
this phase of mass transfer proceeds, we can map out the evolutionary trajectory of various 
sytcms in the log(r/inorm) — log/ diagram if we assume that the system's total mass is 
conserved and the donor's radius is marginally in contact with its Roche lobe. By way 
of illustration, the bottom panel of Figure 1 shows two stable, conservative mass-transfer 
(GMT) trajectories: The (blue) dashed trajectory is for a system of mass Mtot = 1.4 Mq; the 
(pink) dotted trajectory is for a system of mass Mtot = 0.8 Mq. We have assumed that both 
of these systems began the mass-transfer phase of their evolution with an initial mass ratio 
go = 2/3. Hence, the starting point of both trajectories lies on the termination boundary for 
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inspir ailing systems having mass ratios of = 2/3. For systems with Mtot > Met, the mass 
of the accretor will exceed M^h when q drops below the value, 

5,h = ^-l, for Mtot >Meh. (14) 

With the expectation that something catastrophic {e.g., a Type la supernova explosion) will 
occur when this happens, it is reasonable to assume that CMT trajectories with Mtot > M:h 
will terminate at a point in the amplitude- frequency diagram that is marked by ^ch- The 
locus of points that is defined by the termination points of these trajectories defines another 
interesting astrophysical boundary in LISA's "absolute" amplitude-frequency diagram. This 
termination boundary has been drawn as a thick, (green) dashed curve in the bottom panel 
of Figure 1. 

The inspiral trajectory drawn for K — 0.813 (Mtot — 2.4 Mq) and the curve marking 
the termination of various inspiral trajectories in the top panel of Figure 1 define boundaries 
in the amplitude-frequency domain outside of which no DWD system should exist if it has a 
mass ratio g < 2/3. Analogous domain boundaries can be constructed readily for other values 
of q. For each value of q, the shapes of the bounding curves are roughly the same as shown 
in the top panel of Figure 1, but for higher (lower) values of q the right-hand termination 
boundary shifts to higher (lower) frequencies and the limiting inspiral trajectory (set by the 
parameter K) shifts to higher (lower) strain amplitudes. Given our present understanding of 
the structure of white dwarfs, it seems extremely unlikely that any DWD binary systems can 
exist outside of the domain that is defined by the bounding curves for systems with q — 1 
(see, for example. Figure 3). 




Fig. 1. — Top panel: Example evolutionary trajectories (lines with arrows pointing to the 
upper-right) for detached, DWD systems that are undergoing a GR-driven inspiral are dis- 
played in this log(r/inorm) — log / diagram, where r/inorm is specified in meters and / is spec- 
ified in Hz. The three trajectories represent systems having dimensionless mass parameters 
K = 0.813, 0.474, and 0.271 as indicated; assuming a mass ratio q = 2/3 for all three systems, 
this corresponds to total system masses of 2.4, 1.4, and O.SMq, respectively. The termination 
point of each detached, inspiral trajectory is marked by a plus symbol. The solid red curve 
passing through the plus symbols forms a termination boundary for inspiraling systems that 
have a mass ratio g = 2/3 but a wide range of total masses, O.O6M0 < Mtot < 2.4M0; for this 
mass ratio, 2.4M0 represents the maximum allowed total system mass because Ma = Mch- 
Dotted black lines having a slope of —2 depict various "chirp isochrones;" as indicated, the 
isochones identify systems whose characteristic evolutionary time scales, Tchirp, range from 
10^ — lO^yrs. Bottom panel: In the same log(r/inorm)— log/ parameter space, evolutionary tra- 
jectories (curves with arrows pointing to the lower-left) are displayed for DWD systems that 
are undergoing a stable phase of GMT. The two illustrated trajectories have been drawn for 
systems that begin their mass-transfer evolutions with a system mass ratio go = 2/3 - hence 
the trajectories begin at the termination points of the respective inspiral trajectories shown 
above - but they have different total system masses: Mtot = 1-4M0 (blue dotted curve) 
and Mtot = O.8M0 (pink dotted curve). The (green) thick dashed curve represents a low- 
frequency boundary for semi-detached, mass-transferring DWD binaries because, in order for 
a system to move beyond this boundary, the mass of the accretor M^ would have to exceed 
Mch- The single point marked "A" in both panels [ log(r/inorm), log / ] = [—1.65, —2.15] can 
represent DWD systems having multiple (g, Mtot) pairs. 
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3.2. Time-Dependence 



Up to this point, we have described key features of DWD evolutionary trajectories 
in LISA'S "absolute" amplitude-frequency diagram without referring to the rate at which 
the evolution of any given system proceeds. Here we investigate the time scales on which 
significant changes in various system parameters and, as a consequence, the rate at which 
measurable changes in the GW signature occur. Using Eqs.(7) and (8), and assuming Mtot 
as constant, we can write the time-rate-of-change of the amplitude and frequency, for both 
the inspiral and CMT phase, as follows: 



din /in 



dt 



dt 



dt 



din/ 
dt 



^ ainJorb ^ ^d\nQ 



dt 



dt 



(15) 



During the inspiral phase of DWD binary evolution dlnQ/dt = 0, so the evolution is 
driven entirely by the loss of angular momentum due to gravitational radiation. According 
to Peters & Mathews (1963) (see also Misner et al. (1973)), the time-dependent behavior of 
Jorb is described by the relation. 



-'orb(^) = -^0(1 - t/TchirpY^ 



where the inspiral evolutionary time scale is. 



chirp 



Hence, 



256 G^l, 



d In Jorb 



5 



647r2 \rhP 



dt 



8t, 



(16) 



(17) 



(18) 



chirp 



For the CMT phase, dlnQ/dt ^ 0. Hence both the terms on the right-hand side of 
Eq. (15) affect the evolution. From the work of Webbink & Iben (1987) and Marsh et al. 
(2004), we deduce that during a phase of stable CMT (see Kopparapu (2006) for a detailed 
derivation) , 

dt ~ UAcJrebirp' ^''^ 
where A^(Mtot, q) = Cd ~Crl ^ parameter that is of order unity. The quantities (d and (^.i^ 
represent the change in the donor's radius and Roche lobe radius, respectively, as a function 
of its mass (Marsh et al. 2004)^. It should be emphasized that the timescale on which DWD 



''Here we use the notation to represent the change in the radius of the donor with respect to its mass, 
whereas Marsh et al. (2004) use (2 indicating the donor as the secondary star. Also, they use a slightly 
different definition for The relation between their's (Crt) and our's (Crl) is : CrL = Crl + 2(1 — q)- 
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systems evolve during both the inspiral and CMT phases is ~ Tchirp, as indicated by Eqs. (18) 
and (19). It is for this reason that we have drawn various "chirp isochrones" in both panels 
of Figure 1; as can be ascertained from Eq.(17), each isochrone depends only on the product 
of rh and and, hence, has a slope of —2 in the figure panels. 

In practice, for a given source, LISA will be unable to measure changes in the strain 
amplitude h at the levels predicted by expression (15) because variations in h do not ac- 
cumulate secularly over time and, in particular, will not contribute to the observed phase 
of the signal. Since, as indicated in Eq.(l), the phase depends on the time-rate-of-change 
of the GW frequency we will concentrate on the dlnf/dt expression in Eq. (15) from here 
onwards. Combining Eqs. (15), (18) and (19), a concise form of this expression can be 
written to indicate its behavior during both the inspiral (Nelemans et al. 2001) and CMT 
(Nelemans et al. 2004) phases as: 



df 3/ 



dt 8tc 



chirp 



2g-l 



(20) 



where. 



g — (inspiral phase); (21) 

9 = (CMT phase). (22) 

Since depends on both Mtot and q, there exists a critical q = 5crit(M;ot) for which 
AC(Mtot, q) — 0. For systems with q > q'crit(M;ot); a phase of unstable mass transfer ensues 
and our present analysis becomes invalid in that regime. Hence gcrit(M;ot) represents the 
hmiting mass ratio for a system to be in stable CMT phase. 



4. Detectability of DWD Systems 
4.1. Systems with Non- negligible Frequency Variations 



As we have discussed, the physical processes that drive the evolution of DWD binaries 
operate on a "chirp" timescale, and Tchirp is typically much longer than the operational 
time for LISA (assumed one year here). Hence, the time- variation of a given system's GW 
frequency f(t) can be well approximated by a truncated Taylor series expansion in time and, 
using Eq. (1), the observed phase of the GW signal 4>o can be written in the form 



^o{t) = 00 + 27r/ot + 27r 



(23) 
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where /o is the signal frequency at time t — 0, / = df/dt, f — d'^f/dt^, and so on. If we 
truncate the Taylor series at / and this observed signal (O) is compared to a computed 
template (C) that assumes a continuous- wave signal and therefore has a phase that increases 
only linearly with time, 0c (^) = (0o + 27r/ot), the amount of time for the 0-C phase difference 
to reach ■n/2 will be, 

to-c = (2 l/l)-^/^ (24) 
Substituting for / from Eq. (20), we obtain 



to- 



V3|l-2alfn/ 



1/2 



2^|/o 



487r2|l-2^|Vr/i, 



1/2 



(25) 



If the function g in Eq.(25) is independent of h and / — as is the case for the inspi- 
ral phase of DWD evolutions — then curves of constant to-c in the amplitude-frequency 
diagram will be straight lines having a slope of —3. An example of this curve, assuming 
to-c = 1 year, is shown in Figure 3 joining the two low frequency points. Any inspiral 
system that lies to the right of this constant to-c line will lose phase coherence in less than 
one year of observation if one assumes that it emits continous-wave radiation. An analogous 
one-year demarcation boundary can be drawn for DWD binaries that are undergoing a phase 
of stable, CMT by evaluating Eq. (25) using the function g{q,M^^ot)- Because this function 
generally is of order unity, the one- year demarcation boundary for mass-transferring systems 
is generally well-approximated by the line segment that marks the one-year demarcation 
boundary for inspiral systems. LISA will be able to measure frequency evolution for DWD 
systems that lie to the right of the to-c = 1 year line and, as will be discussed in the follow- 
ing section, it will then be possible to determine the distances and masses for these systems. 
The downside is that millions of DWD systems will lie to the left of the one year to-c curve 
in the low frequency region, for which frequency evolution cannot be measured. 



4.2. Determination of Distances and Masses 

An analysis of a one- year-long LISA data stream that utilizes a proper set of frequency- 
varying strain templates should be able to determine the rate at which the GW frequency 
is changing in both inspiral and mass-transfering DWD systems. When used in conjunction 
with the measurement of /inorm and /, an accurate measurement of / for any source should 
permit a determination of the distance to the source r and should give information about 
the chirp mass Ai or the individual component masses of the binary system, as follows. 

Equation (10) provides a relation between the three unknown binary system parameters 
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r, Mtot and g, and the experimentally measurable parameters / and /inorm, namely, 



ML 



L(l + g)2j 



-.12 



norm 



r 



(26) 



A second relation between the unknown astrophysical parameters and measurable ones is 
provided by combining the expression for / in Eq. (20) with the definition of TcMrp given in 
Eq. (17). Specifically, we obtain. 



r(l - 2g) 



5c 
24^ 



/ 



"'normj 



(27) 



With only two equations, of course, it is not possible to uniquely determine all three of the 
binary's primary system parameters. However, in the inspiral phase = 0, so Al and r can 
be determined as was shown in Schutz (1986). 

During the CMT phase of an evolution, the function g{M^^i,q) is nonzero so Eq. (27) 
does not provide an explicit determination of r. However, the requirement that Rd — Rl 
provides an important additional physical relationship between the unknown astrophysical 
parameters and measurable ones. Specifically, by setting R^ from Eq. (12) equal to R^ from 
Eq. (13) and using Kepler's law to write a in terms of /, we obtain. 



R% 



GM, 







1/2 



/ 



7r2(0.0114)' 



ch 



Ma 







-1/2 M, 



tot 



M, 







(^)i7(M,,g), (28) 



where. 



H{Ma,q) = ( 



1/2 



0.49 g2/3 



0.6 g2/3 + in(i + gl/3) 



3/2 



1 - 



M^\4/3 



-3/4 



X 



1 + 3.5 



-2/3 



+ 



-In 



(29) 



Hence, taken together, Eqs. (26)- (28) can be used to determine all three primary system 
parameters - r, Mtot, and q - from the three measured quantities, /inorm, /, and /. 
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-2.5 
Log[f(Hz): 



Fig. 2.— Plot of logT versus log/, where T = {[-ff/f^^V^^^, for stable CMT trajectories 
having five different system masses ranging from 0.6Mq to 2.0Mq, as labeled. Two trajec- 
tories are shown for each value of Mtot- Dashed curve is the analytic solution provided by 
Eqs. (30) and (32); solid curve shows the numerical solution obtained from Eqs. (26) through 
(28). Different values of q are identified by various symbols along each trajectory; note that 
the trajectories with Mtot = 2M0 and I.SMq terminate at values of g = Qch = 0.38 and 0.25, 
respectively. 



- 15 - 



We are unable to solve this set of equations analytically due to the complexity of the 

functions g{Mtot,(l) and H{Md,q). However, the formulae that Paczynski (1967) adopted 
for R(i{M(i) and Rl{q) lead to much simpler expressions for both of these functions, namely, 
g = [|(1 — g)/(2 — 3q)] and H = 1. In this case Eqs. (26)-(28) can be combined to give q in 



terms of / and /: 



g2(l + g)(l-^g)' 



2123 fw 



53c 



15 



(30) 



where a = 0.0141 (GM0i?|)^/2 Qnce q is known, r can be obtained using Eq.(27) in con- 
junction with Paczyhski's g{q) relation; then Mtot can be obtained from Eq.(26). 



5c 



-/ 



Ml 



247r2 Ih 
5^c 



tot 



norm 
3U5 



P 



(2-3g); 



L215 . a-VSGSJ 



1/5 



/ 



11 



(31) 
(32) 



For any Mtot < 2Mch, these three equations are valid for mass ratios over the range < q < 
2/3 because, for Paczyhski's model, ^crit = 2/3 independent of Mtot- 

The solid curves in Figure 2 illustrate results obtained numerically from a self-consistent 
solution of Eqs. (26)- (28); the dashed curves illustrate results obtained analytically from 
expressions (30) and (32). Across the parameter domain defined by the two observables 
log(/) and log(r) — where 



(33) 



is measured in seconds — each curve traces a constant Mtot "trajectory" with q varying 
along each curve, as indicated. At high frequencies, each curve begins at a value of q that 
is slightly below ^crit; a-t low frequencies, the curves have been extended down to q — 0.05, 
unless Mtot > M^h, in which case the curve has been terminated at the value q — qch, as 
given by Eq. (14). The general behavior of these curves can best be understood by analyzing 
analytic expression (30). Over the relevant range of mass ratios < g < 2/3, the analytic 
function. 



anal 



0.0521 



3 ^ 3- 
q'{l + q){l--q 



-1/10 



seconds , 



reaches a minimum value (Fmin = 0.077 seconds) when q — ^extreme; where 



^extreme 



0.2836. 



(34) 



(35) 



Moving from high frequency to low frequency along each Mtot "trajectory," the function F 
steadily drops until q — ^'extreme and F = Fmin- (This behavior holds for the solid curves 
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as well as the dashed curves, although the precise values of Fmin and ^extreme are different 
for each solid curve.) When q drops below ^extreme [based on the function Qch, this will only 
happen along curves for which Mtot < (1 + g'extreme)^ch = l-SSM©], each curve climbs back 
above Fmin, reflecting the fact that Eq. (30) admits two solutions over the relevant range 
of mass ratios. This, in turn, implies that for mass-transferring DWD systems that have 
log(/) < —1.74, a measurement of / will generate two possible solutions - rather than a 
unique solution - for the pair of key physical parameters (Mtot, ?)• 

Once LISA has measured / as well as / for a given DWD system. Figure 2 provides 
a graphical means of determining the values of Mtot and q for the system, assuming it 
is undergoing a phase of stable GMT. We do not expect that LISA will probe the entire 
parameter space depicted in this figure, however. As discussed above, we expect that LISA 
will only be able to detect frequency changes in systems for which io-c ^ 1 yr. Using 
expression (24), this means that LISA will only be able to measure / for systems that have, 

F > 2.57 X IQ-^f-^/^ seconds . (36) 

The dashed black fine in Figure 2 with a slope of —8/5 that is labeled "io-c = 1 year" shows 
this boundary; the parameter regime that can be effectively probed by LISA lies above and 
to the right of this line. 



5. Conclusions 

5.1. Principal Findings 

Once the distance has been determined to individual LISA sources, as was outlined in 
§4.2, it will be possible to place them in an "absolute" amplitude-frequency diagram, that is, 
in a log(r/inorm) — log(/) diagram. The location of individual sources in such a diagram should 
help us understand a great deal about our Galaxy's DWD population. Our expectation is 
that systems in different evolutionary states will fall into several distinct sub-domains of 
LISA'S "absolute" amplitude-frequency diagram, and the diagram will exhibit natural zones 
of avoidance as well. Our models of inspiral and GMT systems permit us to predict where 
the boundaries will lie between these various population sub-domains. As depicted in Figure 
3, the principal sub-domains can be identified as follows. 
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Fig. 3. — DWD population boundaries are drawn in LISA's "absolute" amplitude-frequency 
diagram; five boundary curves (identified as curves A-E in the text) divide the diagram 
into four principal population sub-domains, numbered I-IV. Only inspiraling systems will 
be found in Region I; systems in a phase of stable GMT will only be found in Region II; 
Region III contains inspiraling systems that will eventually encounter a phase of unstable 
mass transfer; and systems in Region IV are progenitors of Type la supernovae. Zones of 
Avoidance: As indicated, no DWD systems will be found in the regions that lie outside of 
the two thick red curves (identified in the text as boundary curves A and B). The straight 
line of slope —3 drawn in the lower left-hand quadrant (identified as curve T in the text) is 
the to-c = 1 yr demarcation line; LISA will only be able to measure / values for systems 
that lie above and to the right of this line. 
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• Zones of Avoidance: No DWD systems will be found in the region above the boundary 
line that is defined by an inspiral trajectory with K — 1; according to Eq. (10), this 
boundary is defined by the expression, 

log(r/inorm) = 0.731 + ^ log / . Boundary A (37) 

Also, no DWD systems will be found in the region below the bounding curve that is 
defined by the inspiral termination boundary for systems with q — 1; over the region 
of parameter space shown in Figure 3 this bounding curve is given approximately by 
the expression, 

log(r/i„orm) ~ 0.703 + 0.637 log/ -0.017 (log /)2 

+ 0.298 (log /)^ + 0.061 (log /)^ Boundary B (38) 

Both of these population boundaries have been drawn as solid red curves in Figure 3. 

Detached, inspiraling DWD systems will be found throughout the domain whose upper 
and lower borders are defined, respectively, by boundary curves A and B. This region of 
parameter space can be subdivided into the following two principal population sub-domains. 

• Region I: Only inspiraling DWD systems will be found in the region of parameter space 
that is bounded above by curve A and below by the locus of points that identifies semi- 
detached systems for which Ma = M^h (the green solid curve in Figure 3); over the 
region of parameter space shown in Figure 3, this curve is given approximately by the 
expression, 

log(r/inorm) ~ 0.761 + 1.005 log / + 0.700(log + 0.700(log ff 

+ 0.214(log/)^ + 0.023(log/)^ Boundary C (39) 

• Region II: Mass-transferring DWD systems will only be found in the region of param- 
eter space that is bounded above by curve C and below by curve B. 

Region II can be further subdivided into two significant population sub-domains as follows: 

• Region III: An inspiraling DWD system that has q > ?crit(M;ot) value will encounter 
an unstable - rather than a stable - phase of mass-transfer when the less massive star 
initially makes contact with its Roche lobe. The region of parameter space in which 
these systems will be found at the onset of mass-transfer is bounded below by curve B 
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and above by two intersecting curve segments: At log / > —1.38, curve C defines the 
upper boundary; for log / < —1.38, the relevant upper boundary has been drawn as a 
blue solid curve in Figure 3 and is given approximately by the expression, 

log(r/inorm) ~ 2.141 + 1.686(log/) 

- 0.141(log/)^ + 0.007(log/)^ Boundary D (40) 

• Region IV: DWD systems that enter a phase of stable CMT with a total system mass 
Mtot > Mch will be found in a region of parameter space that is bounded above by 
curve C and below by two intersecting curve segments: For —1.55 < log / < —1.38, 
curve D defines the lower boundary; for log / < —1.55, the relevant lower boundary 
has been drawn as a pink solid curve in Figure 3 and is given approximately by the 
expression, 

log(r/inorm) ~ " 1.381 - 2. 108 (log /) 

- 1.394(log/)^-0.167(log/)^ Boundary E (41) 

Boundary curve E is defined by the stable CMT trajectory for a system with M^ot = 

An additional demarcation line has been drawn in Figure 3 that is associated with the 
to-c = 1 yr boundary. This "time boundary" line T is defined by setting g — and ^o-c = 1 
yr in Eq. (25), that is, it is given by the expression, 

log(r/i„orm) = -8.498 - 3 log / . Boundary T (42) 

LISA will be unable to determine distances to DWD systems that lie below and to the left 
of this demarcation line because their orbital periods and associated GW frequencies are 
changing so slowly that LISA will be unable to measure with confidence the value of / for 
these systems. Studies of DWD populations will therefore benefit most from the data that 
LISA collects on systems that lie above and to the right of boundary fine T. The segment 
of this line that bounds Region II has been drawn as a dashed line to emphasize that it is 
only an approximate one year boundary for mass-transferring systems. 



5.2. Discussion 

The particular boundaries A-E of the population domains that are depicted in Figure 3 
arise as a consequence of the specific mass-radius relationship (Eq. 12) and, hence, the equa- 
tion of state that we have chosen to use to describe the structure of individual white dwarfs. 
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Because it is generally believed that the mass-radius relationship given in Eq.(12) represents 
the properties of white dwarfs quite well, we expect that the population sub-domains and 
zones of avoidance shown in Figure 3 will map well onto LlSA's observationally determined 
"absolute" amplitude-frequency diagram. Systems that are found on the "wrong" side of 
a given boundary - for example, any system that lies in one of the zones of avoidance, or 
CMT systems that are found in Region I - will be of particular interest because they may 
provide evidence that the equation of state that we adopted is not general enough to properly 
describe this stellar population Two examples suffice to illustrate this point. First, in the 
presence of tidal stresses, the donor stars in mass-transferring DWD systems may be hotter 
than assumed here (Bildsten 2002; Deloye & Bildsten 2003), which will affect the frequency 
at which the donor fills the Roche lobe. Consequently, the inspiral termination boundary 
B in Figure 3 will change. Second, because individual neutron stars obey an entirely differ- 
ent mass-radius relationship and generally seem to have masses > M^h, inspiraling double 
neutron-star systems will likely be distinguishable from DWD systems because they will lie 
in the zone of avoidance above boundary curve A, as defined above. 

DWD systems that are undergoing a phase of stable, CMT and that are found to reside 
in Region IV of the "absolute" amplitude-frequency diagram can be identified as progenitors 
of Type la supernovae. Efforts to better understand the origin of supernova explosions in 
old stellar populations will especially benefit from follow-up studies that identify the optical 
(or UV or x-ray) counterparts to these Region IV systems. Inspiraling systems that are 
found to reside in Region III may prove to be equally interesting candidates for follow-up 
studies. Detached DWD systems in Region III are destined to enter a phase of unstable 
mass transfer (likely accompanied by super-Eddington accretion, see Gokhale et al. (2006)) 
that will significantly transform the system's properties on a dynamical, rather than a chirp, 
timcscale. Such rapid mass-transfer events may lead to merger of the binary components, 
perhaps followed by an explosion. 

As stated above, we are confident that the population sub-domains and zones of avoid- 
ance shown in Figure 3 will map well onto LlSA's observationally determined "absolute" 
amplitude-frequency diagram. We are less confident about the degree to which LlSA's mea- 
surements of / will match the values that are predicted by our simplified model of the 
slow, orbit-averaged time-evolution of DWD systems (as displayed, for example, in Figure 
2). Mass-transferring binary systems, in particular, are notoriously messy laboratories. For 
example, significant and unexplained variations in the mass-transfer rate can arise in an 
individual system over times that are much shorter than the GR-driven evolutionary time 
scale; this can introduce significant short-term variations in /. Also, magnetic fields can be 
effective at carrying away mass and angular momentum from a system, thereby violating 
our assumption of conservative mass transfer. V407 Vul and RX J0806-I-1527 (Marsh & 
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Nelemans 2005) provide perhaps the best examples of the type of confusion that is hkely 
to arise when attempts are made to extract measurements of / from a one-year-long LISA 
data stream. These are optically identified, mass-transferring binaries that are thought to 
be DWD systems because their orbital periods are less than ten minutes. The best available 
measurements of period variation in both of these systems indicate that their orbits are 
slowly shrinking, rather than slowing growing larger as would be predicted by our model of 
stable CMT. Although, in the mean, the long-term evolutionary behavior of these systems 
is likely to agree with the predictions of our model, fluctuations about the mean that occur 
on a timescale that is short compared to Tchirp may totally confound our ability to interpret 
LISA'S measurements of /. On the bright side, if a significant number of LISA sources ex- 
hibit noticeable deviations away from the mean behavior predicted by our simplified model, 
in the end we are likely to gain a more complete understanding of the evolution of such 
systems. 

Wc acknowledge numerous very useful discussions that we have had with J. Frank, V. 
Gokhalc, P. Motl, and X. Peng about the evolutionary behavior of DWD systems. We 
appreciate the encouragement that we received from S. Larson and M. Benacquista to push 
this work to completion after they reviewed a preliminary presentation of our population 
boundary diagram. We also appreciate the advice of an anonymous referee, who urged us 
to shorten the original manuscript making our key point more succinctly. This work has 
been supported, in part, by NASA grant NAG5-13430 and by NSF grants AST-0407070 and 
PHY-0326311. 
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